Cyclotron radiation and emission in graphene 
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Peculiarity in the cyclotron radiation and emission in graphene is theoretically examined in terms 
of the optical conductivity and relaxation rates to propose that graphene in magnetic fields can be 
a candidate to realize the Landau level laser, proposed decades ago [H. Aoki, Appl. Phys. Lett. 48, 
559 (f986)]. 



Introduction — There has been an increasing fasci- 
nation with the physics of graphene, a monolayer of 
graphite, as kicked off by the experimental discovery of 
an anomalous quantum Hall effect(QHE).[l, 2] The fasci- 
nation comes from a condensed-matter realization of the 
massless Dirac-particle dispersion at low energy scales 
on the honeycomb lattice, [l|, [3, 0, S, which is behind 
all the peculiar properties of graphene. In magnetic fields 
this appears as unusual Landau levels, where (i) the Lan- 
dau levels (= y/nhojc, n: Landau index) are unevenly 
spaced, (ii) the cyclotron frequency uJc = {2e/ cTt^vfVB 
is proportional to ^/B rather than to B, and (iii) there is 
an extra Landau level right at the massless Dirac point 
[E = 0), which is outside the Onsager's semiclassical 
quantization. 0] While various transport measurements, 
as exemplified by the quantum Hall effect, have been ex- 
tensively done, optical properties are also measured. For 
example, Sadowski et al have performed a Landau level 
spectroscopy for a large graphene sample. Inter-Landau 
level transitions are observed at multiple energies, which 
is due to a peculiar optical selection rule (|n| <-> |n| -|- 1 as 
opposed to the usual n <-> n + 1) as well as to the uneven 
Landau levels. 

Now, if we look at the QHE physics, cyclotron emis- 
sion from the QHE system in non- equilibrium has been 
one important phenomenon. Experimentally, this typ- 
ically appears as a strong cyclotron emission from the 
"hot spot" , a singular point in a Hall-bar sample where 
the convergence of electric lines of force puts the elec- 
trons out of equilibrium, [sj Theoretically, one of the 
present authors proposed a "Landau-level laser" for non- 
equilibrium QHE systems. [9j The basic idea is simple 
enough: we can exploit the unusual coalescence of the 
energy spectrum into a series of line spectrum (Landau 
levels) realize a laser from a spontaneous emission if we 
can make a population inversion, where the photon en- 
ergy (= cyclotron energy in this case) is tunable and falls 
on the terahertz region for B ^ lOT. However, the most 
difficult part is the population inversion, since if we e.g. 
optically pump the system, the excitation would go up 
the ladder of equidistant Landau levels indefinitely. 

This has motivated us to put a question (Fig[T]): will 
the graphene Landau levels, with uneven spacing among 
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FIG. 1: (Color online) Cyclotron absorption (green) and emis- 
sion (red) processes schematically depicted for the ordinary(a) 
and graphene(b) quantum Hall systems. Black lines represent 
the band dispersion, while blue lines Landau levels. 



their peculiarities, favor in realizing such a population 
inversion? In this Letter we show that this is indeed the 
case, by actually calculating the optical conductivity as 
well as the relaxation processes. The message here is 
graphene is a candidate for the Landau-level laser. 

Optical conductivity in graphene — Low-energy 
physics around the Fermi energy in graphene is described 
by the massless Dirac Hamiltonian,[4] 
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where fi? is the velocity at Ep, tt* = tt^; ± jtTj,, tt = 
p -f eA, A the vector potential representing a uniform 
magnetic field B = rot A, and the 4x4 matrix is spanned 
by the chirality and (K, K') Fermi points. In magnetic 
fields the energy spectrum is quantized into Landau lev- 
els, 
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for a clean system, where n = 0, ±1, 
index. 



is the Landau 
^yti/eB the magnetic length. Here we 
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consider realistic systems having disorder with the self- 
consistent Born approximation (SCBA) introduced by 
Andodli^ to calculate the optical conductivity. 
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FIG. 2: A typical density of states for graphene for an inter- 
mediate disorder (F/ojc = 0.03 here). 



The optical conductivity is given by 
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X [Tr (j„IniG(e)i0(G+(e + hcu) - G+{s))) 
-Tr (j„(G-(e) - G-(e - hu;))jplmG{e))] 



(4) 



where a, (3 = x, y, f{e) the Fermi distribution, and G^ = 
G(e ± «(5). For Green's function G, with the self-energy 

E„(e) =rE„,[£ -sgn(n)y^ in the SCBA, 
the Landau level broadening is given by F = tioVq if 
we assume for simplicity short-ranged random potential, 
V = J2i^oS{'r — Ti). The light absorption rate is then 
related to the imaginary part of the dielectric function, 
e{u}) = 1-1- iaxxi^^)/£o^, so we can look at Keaxx{to)- 

In order to discuss the optical conductivity in graphene 
we need the current matrix elements across Landau lev- 
els. The eigenfunctions of the Hamiltonian ([1]) dictate 
an unusual selection rule, |?tJ_- \n'\ = ±1 in place of the 
ordinary n — n' = ±1, with[10| 

Jx'" = vpCnCn' [sgn(n)(5|„|_ij„/| + sgn(n')(5|„|+i,|„/|] , 
= ivpCnCn' [sgn(n)(5|„|_i_|„,| - sgn(n')(5|„|+ij„'|] ,(5) 

where G„ = l(n = 0) or l/\/2 (otherwise). 

We have numerically obtained the Green's function 
and optical conductivity. While in usual cases the broad- 
ened Landau levels are uniformly merged or separated as 
F is varied, there is a striking difference for graphene, 
where the Landau levels (oc ^/n) are unevenly spaced, so 
that the broadened Landau levels overlap to lesser ex- 
tent as we go to the central one {n 0), as typically 
depicted in Fig[2l Namely, for an intermeditate value of 
T/uJc only the n — Landau level stands alone while the 
other levels form a continuous spectrum. 

We now look at the optical conductivity in Fig[3]for the 
Fermi energy at ep = (energy for the Dirac point), each 




FIG. 3: (Color online) Optical conductivity against uj is shown 
for temperatures kBT/TtWc = — 2.0, for a fixed value of 
V/TiLOc = 0.01 with = 0. 



resonance peak can be assigned to an allowed transition 
with the selection rule (eqn®). The largest peak around 
lo/lOc = 1 corresponds to the transition between n = 
±1, while the peaks at higher frequencies come from the 
transition across the Fermi energy, —n ^ n ± 1. If we 
turn to the temperature dependence in the figure, we 
immediately notice a peculiar phenomenon: there is a 
peak, in the region lj/lOc < 1, that grows, rather than 
decays, for higher T. We can identify this as coming from 
the unusual Landau levels in graphene: As T is raised 
with the Fermi distribution function becoming longer- 
tailed, higher Landau levels begin to be occupied, which 
enables the transitions among higher Landau levels, n ^ 
nil, to take place. While this would not cause new 
lines to appear for equidistant Landau levels, this does 
so for the unequally spaced Landau levels (cx \/|n[) for 
uj/loc < 1 transitions. So we can identify this property 
as one hallmark of the "massless Dirac" dispersion. 

Previously, the ogtical conductivity has been obtained 
by Gusynin et al, llj I4I who have derived the analyti- 
cal expression for the optical conductivity, but the self- 
energy from the disorder was set to a constant, while 
we have calculated the self-energy self-consistently with 
SCBA. Sadowski et al.[13[ also presented a similar ex- 
pression for the conductivity with a constant self-energy 
as well. The present result qualitatively agrees with 
these, but the new findings here are, first, the full de- 
pendence on the kBT/fiiOc, including the growing of low- 
frequency peaks at low temperatures. Secondly, we point 
out that the situation as depicted in Fig[2] should be in- 
teresting for the cyclotron resonance and emission in non- 
equilibrium situations induced by e.g. an optical pump- 
ing with laser beams. Namely, the electrons excited to 
higher energies will relax down to the n = \ level across 
the continuum spectrum, so that the population inversion 
across n = and n > 1 should be easier to be realized. 
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Relaxation processes — To quantify this idea, we have 
to consider the relaxation processes which should control 
the population inversion. For the ordinary quantum Hall 
systems the relaxation processes have been extensively 
discussed. Specifically, Chaubet et al. 14 1 1^ discussed 
dissipation mechanisms, where spontaneous photon ra- 
diation and coupling with phonons are examined on the 
basis of Fermi's golden rule. Other dissipation processes 
such as electron-electron scattering or impurity scatter- 
ing, which conserve the total energy, do not contribute 
to inter-Landau level processes in the absence of external 
electric fields (while Chaubet et al. have focused on ef- 
fects of finite electric fields in the QHE breakdown where 
inter-Landau level processes are involved). 

So we extend the discussion by Chaubet et al. to re- 
laxation processes in graphene. We first estimate the ef- 
ficiency of the photon emission with Fermi's golden rule: 
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Here \i){ei) is the wavefunction (energy) in the initial 
state while f stands for the final states, and Hint the 
electric dipole interaction between the electromagnetic 
field and electrons. When the wavelength of light is much 
larger than the cyclotron radius, as is usually the case, 
we have 



W, 



2n V 



2eoVLj 



X d(hLO -I- £/ — ei) = 4q; I I 07, (6) 

where c is the velocity of light, a — e^/ (47reo?ic) the fine- 
structure constant, and we put huj — ef — ei to be the 
cyclotron energy hujc- 

A peculiarity of graphene appears in the current ma- 
trix element (cqn®), for which the rate of spontaneous 
emission, with + 1)| = C„C„+iVf for graphene 

plugged in, reads 



Wgraphene 




(n = 0), 
(n^O). 



(7) 



This expression, another key result here, shows that the 
spontaneous emission rate depends linearly on the cy- 
clotron energy and quadratically on the Fermi velocity. 
This is in sharp contrast with the ordinary QHE systems 
such as the two-dimensional electron gas (2DEG) real- 
ized at e.g. GaAs/AlGaAs interfaces. In this case the 
velocity matrix element + 1)1"^ = (n + l)hLu/2m* 

should be plugged in eqn.Q, which yields 
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This reveals a dramatic difference between graphene and 
usual 2DEG, where the emission rate in the latter is pro- 
portional to the square of the cyclotron energy. 



We can quantitatively realize the difference: The cy- 
clotron energies are 



huj 



heB/m* - 
VF^2heB : 



L7meV 
- 37meV 



(GaAs), 
(graphene). 



for B = 1 T, where we have adopted the value of 
graphene Fermi velocity vp = 1-06 x 10^ m/s,[l3] and 
the GaAs effective mass m* ~ 0.067me. Hence the cy- 
clotron energy in graphene is orders of magnitude larger 
since it scales as ^/B reflecting the Dirac dispersion, while 
the energy is usually proportional to B. If we plug these 
in eqns. ()7|8p . we end up with 



cx B2 ~ 6 X ^(s-i) (GaAs), 
cx ^/B ~ 1 X 10^(s^^) (graphene). 



where the second term in each line indicates the B- 
dependence, while the third term a numerical value for 
_B = 1 T. A conspicuous difference, oc B^ in the former 
and cx ^fB, should sharply affect the behavior. Thus the 
spontaneous photon emission rate is orders of magnitude 
enhanced in graphene in moderate magnetic fields (as in 
the above numbers quoted for B = 1 T.) This indicates 
that the present system is indeed favorable for a realiza- 
tion of the envisaged Landau level laser. 

Now, the dissipation process which competes with the 
photon emission is the phonon emission process, which 
has been discussed for the conventional QHE systems, es- 
pecially in the context of the breakdown of the quantum 
Hah effect [ill . The phonon emission rate is also obtained 
from Fermi's golden rule if we replace the electron-light 
interaction with the electron-phonon interaction. If we 
first consider acoustic phonons, the dissipation rate is 
promotional to the extent of the overlap between initial 
and final wavefunctions both in usual and graphene QHE 
systems, which yields a factor e~^^^) with q the phonon 
wavenumber and £ the magnetic length. In usual QHE 
systems the cyclotron energy is ~ 1 meV and the mag- 
netic length i — y^h/eB ~ 30 nm for B =1 T, while the 
acoustic phonon wavenumber is ~ lA"-'^, so that the over- 
lap factor is exponentially small. The situation is similar 
in graphene, since the magnetic length £ = ^JfijeB is the 
same. So the acoustic phonon emission should be negligi- 
ble in graphene as well in weak electric fields. When the 
applied laser electric field is so intense 1 kV/cm) that 
the Landau levels are distorted and the overlap factor 
grows, the phonon emission may begin to compete with 
the photoemission. 

Are there any other factors that distinguish graphene 
from 2DEG's? In this context we can note that Chaubet 
et al. have further pointed out the following. In an elec- 
tron system confined to 2D a wavefunction has a finite 
tail in the direction normal to the plane, and the phonon 
emission is enhanced through the coupling of the tail of 
the wavefunction and perpendicular phonon modes which 
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propagate normal to the 2D system in the substrate 15 1. 
This way the phonon emission can compete with the 
spontaneous emission in usual QHE systems. By con- 
trast, a graphene sheet is an atomic monolayer, and there 
is only a loose coupling with the substrate. We can also 
consider acoustic phonons coupled with impurity scatter- 
ing, which may compensate the momentum transfer q of 
phonons, and hence the overlap factor e"^''^-* To be 
precise, graphene itself should have phonon modes that 
include the out of plane modes, and their effects is an 
interesting future problem. As for optical phonons, their 
energies are known to be higher than 100 meV for wave- 
length g = in graphene [16|, so that optical phonons 
do not contribute to the dissipation for B ~ a few tesla 
with hui ^ 40 meV. Overall, we conclude that the dissipa- 
tion due to acoustic phonons will be small in graphene in 
the weak electric-field regime. When the pumping laser 
intensity is not too strong to invalidate the present treat- 
ment but strong enough for the population inversion, the 
present reasoning should apply, and we can expect effi- 
cient cyclotron emissions from graphene. 

Entirely different, but interesting is the problem of An- 
derson localization arising from disorder. While this is 
out of scope of the present work, we can expect delocal- 
ized states with diverging localization length at the cen- 
ter of each Landau level are present as inferred from the 
QHE observation, whose detail is an interesting future 
problem. The situation should also depend on whether 
the disorder is short-range or long-range, but, in ordinary 
QHE systems, a sum rule guarantees the total intensity 
of the cyclotron resonance intact, [o'l As for the " ripp les" , 
suggested to exist in actual graphene samples 17| , the 
n = Landau level remains sharp (which is topologically 
protected since the slowly varying potential does not de- 
stroy the chiral svmmetrv[l8l| ) . while other levels become 
broadened, [3| and this favors the situation proposed in 
the presented paper. [2^ 

Summary — To summarize, we have discussed the ra- 
diation from graphene QHE system. We conclude that 
unusual uneven Landau levels, unusual cyclotron energy, 
unusual transition selection rules all work favorably for a 
population inversion envisaged for the Landau level laser. 
An estimate of the photon emission rate shows that the 
emission rate is orders of magnitude more efficient than in 
the ordinary QHE system, while the competing phonon 
emission rate is not too large to mar the photon emis- 
sion. Important future problems include the examination 
of the actual lasing processes including the cavity prop- 



erties, coupling of electrons to the out-of-plane phonon 
modes, etc. We wish to thank Andre Geim for illumi- 
nating discussions. This work has been supported in 
part by Grants-in-Aid for Scientific Research on Prior- 
ity Areas from MEXT, "Physics of new quantum phases 
in superclean materials" (Grant No. 18043007) for YH, 
"Anomalous quantum materials" (No. 16076203) for HA. 



[1 
[2: 

[s: 
[4: 
[s: 

[6: 

[7: 
[s: 

[9: 
[10: 

[11 

[12 
[13 
[14 
[15 
[16 
[17 
[18 
[19 
[20 
[21 

[22; 



K. S. Novoselov et al, Nature, 438, 197 (2005); Nature 
Physics, 2, 177 (2006). 

Y. Zhang, Y. W. Tan, H. L. Stormer, and P. Kim, Nature 
438, 201 (2005). 

J. McClure, Phys. Rev. 104, 666 (1956). 
Y. Zheng and T. Ando, Phys. Rev. B 65, 245420 (2002). 
V. P. Gusynin and S. G. Sharapov, Phys. Rev. Lett. 95, 
146801 (2005). 

N. M. R. Peres, F. Guinea, and A. H. Castro Neto, Phys. 

Rev. B 73, 125411 (2006). 

L. Onsager, Phil. Mag 43, 1006 (1952). 

K. Ikushima, H. Sakuma, S. Komiyama, and K. Hi- 

rakawa, Phys. Rev. Lett. 93, 146804 (2004). 

H. Aoki, Appl. Phys. Lett. 48, 559 (1986). 

T. Ando, J. Phys. Soc. Japan 38, 989 (1975). 

V. P. Gusynin, S. G. Sharapov, and J. P. Carbotte, Phys. 

Rev. Lett. 98, 157402 (2007). 

V. P. Gusynin and S. G. Sharapov, Phys. Rev. B 73, 
245411 (2006). 

M. L. Sadowski, G. Martinez, M. Potemski, C. Berger, 
and W. A. de Heer, Phys. Rev. Lett. 97, 266405 (2006). 
C. Chaubet, A. Raymond, and D. Dur, Phys. Rev. B 52, 
11178 (1995). 

C. Chaubet and F. Geniet, Phys. Rev. B 58, 13015 
(1998). 

S. Piscanec, M. Lazzeri, F. Mauri, A. C. Ferrari, and 

J. Robertson, Phys. Rev. Lett. 93, 185503 (2004). 

J. Meyer, A. Geim, M. Katsnelson, K. Novoselov, 

T. Booth, and S. Roth, Nature 446, 60 (2007). 

Y. Hatsugai, T. Fukui, and H. Aoki, Phys. Rev. B 74, 

205414 (2006). 

A. Giesbers, U. Zeitler, M. Katsnelson, L. Ponomarenko, 
T. Ghulam, and J. Maan, eprint arXiv: 0706.2822 (2007). 
Y. Niimi, H. Kambara, T. Matsui, D. Yoshioka, and 
H. Fukuyama, Phys. Rev. Lett. 97, 236804 (2006). 
Localized states around point defects were detected re- 
cently with STM in B [23|, and their radius was found 
to be comparable with the magnetic length (~ 30 nm). 
We can also note that the n — Landau level in graphene 
has an exactly E — edge mode, of a topological origin, 
right at the center [l^, which may also affect the photon 
absorption / emission. 



